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Since the adoption of the ICH Q8 document concerning the development of pharmaceutical 
processes following a quality by design (QbD) approach, there have been many discussions 
on the opportunity for analytical procedure developments to follow a similar approach. While 
development and optimization of analytical procedure following QbD principles have been 
largely discussed and described, the place of analytical procedure validation in this frame- 
work has not been clarified. This article aims at showing that analytical procedure validation 
is fully integrated into the QbD paradigm and is an essential step in developing analytical 
procedures that are effectively fit for purpose. Adequate statistical methodologies have also 
their role to play: such as design of experiments, statistical modeling, and probabilistic state- 
ments. The outcome of analytical procedure validation is also an analytical procedure design 
space, and from it, control strategy can be set. 
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1. INTRODUCTION 

The concept of quality by design (QbD) has been adopted in the pharmaceutical 
industry through several initiatives such as the FDA's cGMP for the twenty-first century 
(U.S. Food and Drug Administration (FDA), 2007) and process analytical technology (PAT) 
(United States Food and Drug Administration (FDA), 2004) as well as with the regulatory 
documents ICH Q8 (2009), Q9 (2005), and Q10 (2008) and the FDA guidance on process 
validation (U.S. Food and Drug Administration (FDA), 201 1). The general aim is to switch 
from the quality by testing (QbT) paradigm previously implemented in the pharmaceutical 
industry to a development aiming at improving the understanding of the processes and 
products and hence improving product quality, process efficiency, and regulatory flexibility. 

QbD is not new and involves many quality and statistical tools and methods, such 
as statistical designs of experiments, multivariate statistics, statistical quality control, and 
so on. In order to raise the quality of pharmaceutical products, it has been recognized that 
increasing the testing of final products (i.e., QbT) is not adequate (Lionberger et al., 2008). 
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Instead, to increase the quality of pharmaceutical products, quality must be built into the 
products (i.e., QbD) as already done in many other industries. It requires understanding 
how variables involved in formulation and manufacturing processes influence the quality 
of the final product. 

Analytical procedures are also processes, and QbD should also be implemented 
for the development of analytical procedures. Several authors recently stated that QbD 
enables to develop analytical procedures in a systematic and scientific approach (Borman 
et al., 2010; Nethercote et al., 2010; Schweitzer et al., 2010; Ermer, 2011; Rozet et al, 
2013). The understanding and identification of variables affecting method performance is 
achieved at an earlier stage (Borman et al., 2010; Nethercote et al., 2010; Schweitzer et al., 
2010; Ermer, 2011; Rozet et al., 2013). According to ICH Q8(R2) (U.S. Food and Drug 
Administration (FDA), (2011), QbD can be seen as an optimization strategy combining 
design of experiments (DoE) and design space (DS). 

However, developed analytical procedures are not directly usable in laboratories as 
they have to demonstrate that they are indeed fit for their purpose. This demonstration of 
fitness of purpose is generally achieved during the analytical procedure validation phase. 

The aim of this article is to show that analytical procedure validation is fully inte- 
grated into the QbD paradigm and is an essential step in developing analytical procedure 
that are useful for their routine applications. Similar statistical methodologies are also 
implemented in analytical procedure validation among which are DoE and statistical 
modeling. The outcome of analytical procedure validation is also an analytical procedure 
DS, and from it control strategy can be further defined. 

2. ANALYTICAL TARGET PROFILE AND ANALYTICAL PROCEDURE 
VALIDATION 

The development of a QbD compliant analytical procedure starts by the definition 
of its analytical target profile (ATP) that aims at defining the intended purpose of the pro- 
cedure. The ATP compiles a set of characteristics defining what analyte or analytes will 
be measured, in which matrix, over what concentration range(s), as well as the required 
performance criteria of the method together with specifications for these last ones. These 
specifications and characteristics should be linked to the intended purpose of the analytical 
procedure. Examples of ATP for the interested readers can be found in the following refer- 
ences (Borman et al., 2010; Nethercote et al., 2010; Schweitzer et al., 2010; Ermer, 2011; 
Rozet et al., 2012). The requirements of quantitative performances included in the ATP 
are then the validation acceptance limits that must be reached by the analytical procedure 
during the validation phase. In addition, several authors have gone further in the definition 
of the ATP by including maximum acceptable risk of making wrong decisions using the 
results generated by analytical procedures (Rozet et al, 2012). 

3. CRITICAL QUALITY ATTRIBUTE IN ANALYTICAL PROCEDURE 
VALIDATION 

The critical quality attributes (CQAs) of the analytical procedure are the responses 
that are measured to judge the quality of the developed analytical procedures. CQAs are 
defined as "a physical, chemical, biological or microbiological property or characteristic 
that should be within an appropriate limit, range, or distribution to ensure the desired prod- 
uct quality" (ICH Q8, 2009). For chromatographic analytical procedures, the CQAs can 
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be related to the method selectivity, such as the resolution (Rs) or separation (5) criteria 
(Rozet et al., 2013). Other CQAs can be the run time of the analysis, signal to noise ratio, 
the precision and the trueness of the analytical procedure, the lower limit of quantification 
(LOQ), or the dosing range of the analytical method. These CQAs may be directly mod- 
eled through a multivariate (non)linear model. However, in other situations, the modeled 
responses may be different than the CQAs. The CQAs are obtained after the modeling of 
these primary responses. For chromatographic methods, the usual key CQA is resolution of 
the critical pair when optimizing selectivity. However, the resolution depends on the reten- 
tion factor of the two chromatographic peaks involved. Therefore, the retention factors are 
directly modeled instead of the resolution. The resolution can then be computed from these 
modeled responses. 

Nonetheless, CQAs are not limited to separative techniques or only related to the 
qualitative performance of the analytical procedure. Highly important assays in the devel- 
opment and control of pharmaceutical products are quantitative ones. Other examples than 
chromatographic quantitative procedures are immunoassays such as ELISA, q-PCR, rela- 
tive potency assays, and so on. The final aim of any quantitative analytical procedures is to 
provide analytical results of adequate quality in order to make reliable decisions with them. 
Hence, CQAs for quantitative procedures should be at least related to their quantitative 
performances. The validation characteristics that are trueness, precision, linearity, range, 
LOQ of the analytical procedure, and the accuracy of the results obtained by the procedure 
are key CQAs. They should be included into the definition of the ATP together with their 
respective acceptance values. 

The validation phase of any quantitative analytical procedure is therefore fully in 
line with the QbD framework. The CQA that should be monitored during analytical pro- 
cedure validation are measures related to random error (e.g., intermediate precision CV), 
systematic error (e.g., bias or recovery), or the combination of both which is total error. 

4. CRITICAL PROCESS PARAMETERS IN ANALYTICAL PROCEDURE 
VALIDATION 

Analytical procedure validation also involves several factors that are critical pro- 
cess parameters. The first main factor is concentration/ amount/potency range over which 
the procedure is intended to quantify the analyte. This factor is a fixed factor and is 
represented by samples called validation standards or quality control samples of known 
concentration/amount/potency. By opposition, sources of variability that will be encoun- 
tered during the future routine use of the procedure must be included in the validation 
design as random factors, such as operator, equipment, reagent batch, or days. The 
combination of these sources of variability is generally called runs or series. 

5. DESIGN OF EXPERIMENTS IN ANALYTICAL PROCEDURE VALIDATION 

ICH Q8 and FDA guideline highly promote the use of adequate DoE when developing 
pharmaceutical processes. The main designs used in analytical procedure validation are 
nested designs or (fractional) factorial designs or a combination of both. These designs are 
used to estimate the variance components. To have precise estimations, the use of more than 
two levels of each factor is recommended. Nonetheless, the various sources of variation 
included into the analytical procedure validation are generally combined into "series" or 
"runs" to mimic the way analytical procedures are effectively employed routinely. 
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Suppose that for each of the ; concentration level of the validation standards, the 
number of runs is J and that in each run, K replicates are performed. The validation exper- 
iments can then be described, for each of the ; th concentration level studied, by a one-way 
analysis of variance (ANOVA) random model with runs (or series) as random factor: 

X iJk = ^ + a uj + e iJk , a hJ ~ N (0, o£,.) , e ujk ~ N (0, o£ ( ) (1) 

where [ii is the overall mean of the i th concentration level studied of the validation standard, 
iu, + a h j is the mean in run j (j: 1 to J), syj is the residual error, er 2 , is the run-to-run 
variance, and er f 2 ( is the within-run or repeatability variance, both for the ft 1 concentration 
level. 

The overall variability of the analytical method is measured by the intermediate pre- 
cision variance a? p i = er„ ; + a^ t . All these parameters of the variance components model 
can be estimated by restricted maximum likelihood (REML) methods (Searle et al., 1992). 

6. DESIGN SPACE AND ANALYTICAL PROCEDURE VALIDATION 

In the ICH pharmaceutical development guideline Q8 (ICH Q8, 2009), the DS is 
defined as "the multidimensional combination and interaction of input variables (e.g. mate- 
rial attributes) and process parameters that have been demonstrated to provide assurance 
of quality." Therefore, the multidimensional combination and interaction of input variable 
corresponds to a subspace, so-called the DS, where assurance of quality has been proven. 
The main concept lying behind the ICH Q8 definition of DS is assurance of quality (also 
known as quality risk management). It has been already shown that the mean response 
surface obtained during analytical procedure development does not define properly a DS as 
there is no assurance that the CQAs reach their acceptance limits. Instead, probability maps 
answer this DS requirement properly (Peterson and Lief, 2010; Rozet et al., 2013). 

Analytical procedure validation also allows defining a DS: it is the range of concen- 
tration where it has been demonstrated that the procedure provides assurance of quality 
results 

n = P(-k < X - ht < X). (2) 

The objective of the validation phase can be summarized to evaluate whether the reli- 
ability probability it that each future result will fall within predefined acceptance limits (A) 
is greater than or equal to a minimum claimed level 7T m j n (Hubert et al., 2004). The statis- 
tical problem here is twofold: the probability it needs to be estimated, and the uncertainty 
in its estimation must be taken into account when comparing it to :rc m j n . This is not an easy 
problem to solve since it has no exact small sample solution in frequentist statistics. 

Nonetheless, several approaches have been proposed to answer this aim. 

6.1. ^-Expectation Tolerance Intervals 

A first one is to compute ^-expectation tolerance intervals of a defined coverage prob- 
ability (e.g., 95%) at each concentration level of the validation standards using the one-way 
ANOVA random model described in Eq. (1) and comparing it to preset acceptance limits 
as shown in Fig. 1. Using this approach, each future result has at least 95% probability 
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Figure 1 Accuracy profile, depicting at each concentration level of the validation standards the corresponding 
95% P-expectation tolerance intervals (dashed lines). The acceptance limits have been set at ±15% around the 
known concentration values of the validation standards (dotted lines above and below dashed lines). The center 
continuous line shows the relative bias of the assay. The dots are the analytical results of the validation standards 
expressed in relative error values. 



to fall within these acceptance limits. Lebrun et al. (2013) have shown that P-expectation 
tolerance intervals are equivalent to highest posterior density intervals. A non-negligible 
amount of analytical procedures have been validated in such a way (Hubert et al., 2004, 
2007a, 2007b). Figure 1 shows an accuracy profile obtained for the validation of an ana- 
lytical procedure, depicting at each concentration level of the validation standards the 95% 
P-expectation tolerance interval. 



6.2. Out-of-Specification Probability 

Another approach is to estimate the probability to obtain future results outside the 
preset acceptance limits (out-of-specification, OOS). Dewe et al. (2007) has proposed 
to compute this probability for results following the one-way ANOVA random model 
described in Eq. (1). An example is shown in Fig. 2 for the same previous analytical pro- 
cedure. The DS is then the range of concentration over which this probability is smaller 
than a preset maximum value (e.g., 0.05). For each concentration level this probability is 
computed as: 



n *a = p [ X . > ^ _ x j + p [ X . < ; + A j 



= P 



t(f) > 



KRj + \ 
N(R, + 1) _ 



t(f) < 



(li T i + k) - X, 



(3) 



ffi.p.,ijl 



KRi+\ 



N(R,+1) _ 



where / is the number of runs and K the number of replicates by series, N = JK. X t is the 
mean concentration of the results obtained by the method for the i th concentration level and 
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Figure 2 Risk profile, giving at each concentration level of the validation standards the probability to have future 
analytical results falling outside an acceptance value of ±15% around the known concentration values of the 
validation standards, i.e., OOS probability. The maximum OOS probability has been set at 5%. 



<7i p j is the intermediate precision standard deviation for each i concentration level. t(f) 
is a Student distribution with / degrees of freedom computed based on the Satterthwaite 
approximation (Satterthwaite, 1941) and is the ratio between the run-to-run variance and 
the within-run (or repeatability) variance of each concentration level. The use of a Student 
distribution is justified as it is the predictive distribution in this model as demonstrated by 
Lebrun et al. (2013). 



6.3. Continuous Modeling Across Concentration Range: The Bayesian 
Way 

If between-run variances and repeatability variances can be assumed homogenous 
across the concentrations levels of the validation standards, a single linear mixed model can 
be fitted to the validation data, including concentration as a fixed factor. The two previous 
approaches to define a DS can then be extended to this situation. 

A less trivial situation would be to model the analytical procedure results over the 
concentration range in case of heteroscedasticity of between-run and/or repeatability vari- 
ance. Determination of ^-expectation tolerance intervals or probability of OOS in these 
cases could be based on Bayesian approaches as no frequentist solutions are available 
(Rozet etal., 2011). 

In this context model, Eq. (2) is rewritten as the following linear model with random 
slopes and intercepts and residual variance increasing with concentration: 

X ijk = Po+ PlI^TJ + UQj + Uij/lTj + s ijk . (4) 

where the subscripts i stands for the / concentration levels of the validation standards, j 
for the J number of series or runs, and k for the K number of replicates per run. ijlj is 
the I th concentration level of the validation standard and is considered as a reference or 
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Figure 3 Bayesian risk profile, modeling over the concentration range studied the probability to have future 
analytical results falling outside an acceptance value of ±15% around the known concentration values of the 
validation standards, i.e., OOS probability. The maximum OOS probability has been set at 5%. The lower LOQ 
corresponds to the concentration where the OOS probability crosses the maximum OOS probability value of 5%. 



conventional true value. 6= ( „ | are the fixed effects. Additionally, U, = ( M °' J | are the 
random effects of the / runs and are also assumed coming from a normal distribution: 



U; ~iAM0,a B 2 S I. (5) 



Finally, e,^ is the residual error assumed to be independent and coming from a normal dis- 
tribution of variance of. This variance is also given as being dependent on the concentration 
level (. This phenomenon is frequently observed in real-life situations. The general form of 
this variance function is a power of the concentration: 

Oi = a(n T .i) Y - (6) 

Figure 3 illustrates the probability profile for an analytical procedure using this model and 
estimated using Markov Chain Monte Carlo (MCMC) simulations, ft depicts the concen- 
tration range over which the analytical procedure is fit for its purpose. This range represents 
the analytical procedure validation DS. 



7. CONTROL STRATEGY 

QbD development of analytical procedure is useless without defining control strat- 
egy to ensure that the procedure remains under control during its routine application and 
detect deviations. Validation analytical procedure validation also allows defining a control 
strategy using quality control samples. Indeed, the experiments performed allow defining 
^-expectation tolerance intervals that can be used as initial control limits when building 
analytical procedure control charts (Rozet et al., 2007). Out-of-control methods can effi- 
ciently be detected and corrective actions realized by following the daily performances 
of analytical methods on such charts. Indeed, the use of fi -expectation tolerance intervals 
ensures an adequate balance between consumer and producer risks (Rozet et al., 2007). 
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8. CONCLUSION 

Analytical procedure validation fits entirely within the QbD paradigm. In fact, when 
comparing it with pharmaceutical development, analytical procedure validation can be seen 
at the stage 2 of process validation as defined by the recent FDA guideline (U.S. Food 
and Drug Administration (FDA), 2011). Analytical procedure validation is then the perfor- 
mance qualification of the assay. In this context, seeing analytical procedure validation as 
an additional burden in analytical procedure development limited to the "ICHQ2 (ICH Q2, 
2005) check list" exercise should disappear: the validation phase is the confirmation of the 
usefulness of the developed procedure for its future daily application. 
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